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In systems with broken inversion symmetry spin-orbit coupling (SOC) yields a 
Rashba-type spin splitting of electronic states, manifested in a k-dependent split¬ 
ting of the bands. While most research had previously focused on 2D electron 
systems mm, recently a three-dimensional (3D) form of such Rashba-effect was 
found in a series of bismuth tellurohalides BiTeX (X=I, Br, or Cl) |3]. 1^ [5|. [6] . 
Whereas these materials exhibit a very large spin-splitting they lack an impor¬ 
tant property concerning functionalisation, namely the possibility to switch or 
tune the spin texture. This limitation can be overcome in a new class of func¬ 
tional materials displaying Rashba-splitting coupled to ferroelectricity: the fer¬ 
roelectric Rashba semiconductors (FERS) [7l[8]. Using spin- and angle-resolved 
photoemission spectroscopy (SARPES) we show that o-GeTe(lll) forms a prime 
member of this class, displaying a complex spin-texture for the Rashba-split 
surface and bulk bands arising from the intrinsic inversion symmetry breaking 
caused by the ferroelectric polarization of the bulk (FE). Apart from pure sur¬ 
face and bulk states we find surface-bulk resonant states (SBR) whose wavefunc- 
tions entangle the spinors from the bulk and surface contributions. At the Fermi 
level their hybridization results in unconventional spin topologies with cochiral 
helicities and concomitant gap opening. The o-GeTe(lll) surface states (SS) 
and SBR states make the semiconductor surface conducting. At the same time 
our SARPES data confirm that GeTe is a narrow-gap semiconductor with a well 
resolved Dirac point in the 3D-valence band and a giant spin-splitting AkR ^0.11 
of the bulk states. These results show that the o-GeTe(lll) surface/bulk 
electronic states are endowed with spin properties that are theoretically chal¬ 
lenging to anticipate. As the helicity of the spins in Rashba bands is connected 
to the direction of the FE polarization, this work paves the way to all-electric 
non-volatile control of spin-transport properties in semiconductors. 

IV-VI semiconductor compounds such as lead chalcogenides, SuTe or GeTe have been exten¬ 
sively studied by band-structure calculations for which relativistic effects and SOC were in¬ 
cluded ismniin]. Emphasis was primarily put on the band-edge properties and the fundamen¬ 
tal gap m rather than the spin properties. However, recent theoretical work has predicted 
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that in topologically trivial (a-GeTe(lll) the bulk bands are expected to be fully spin-polarized 
with the spin chirality directly coupled to the FE polarization lam- A further development 
of this concept into applications requires first of all an experimental verification and, deeper 
theoretical considerations. 
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Figure 1: Basic properties of (a-GeTe(lll) and experimental geometry, (a) Bril- 
louine zone; (b) RHEED and (c) FEED pattern of the bulk and surface, respectively; (d) 
Te-terminated semi-infinite surface structural model with highlighted distorted rhombohedral 
cell. Orange arrow shows the displacement of the Ge-atom responsible for the polar FE or¬ 
dering. (e) Experimental geometry with schematic drawing of our SARPES experiment, (f) 
Spin-resolved initial state Bloch spectral functions calculated for bulk Brillouine zone along ZU 
and ZA directions, (inset) ~^x,y quiver plot along the green arrow. Calculated spin-integrated 
(g) and spin-resolved (h) BSF from a semi-infinite surface, (i) Logarithmic false color map of 
ISM determinant demonstrating bulk, surface and SBR character of calculated ARPES spec¬ 
tral intensities [TBIIT3] . The green arrow indicates the ky momentum span where bulk Rashba 
splitting is best observed. 

Experiments were performed on 200 nm thick (a-GeTe(lll) monocrystalline films grown by 
molecular beam epitaxy on BaF 2 (lll) substrates [121 [E]? with a Te-terminated surface pre¬ 
pared in situ [see Methods and Supplemental Information (SI) 1,11.1]. Fig. l(b-d) show the 
surface structural model together with results of the sample characterization by in-situ low- 
energy electron diffraction (FEED) and reflection high-energy electron diffraction (RHEED). 
The lack of surface reconstruction in FEED and the presence of Kikuchi lines in RHEED indi¬ 
cate superior quality (a-GeTe(lll) surface with spontaneous FE moments [13]. Below Tc = 700 
K there is a spontaneous paraelectric to FE phase transition in which the cubic rock-salt phase 
changes to a polar rhombohedral structure, forming a FE-phase induced by elongation of the 
unit cell along the (111) direction (see the orange arrow in Fig. Id). A concomitant shift in the 
Te-sublattice is inducing an almost perfect collective FE-order, up to ^92% oriented along the 
(111) direction perpendicular to the sample surface (SI H.2). The resulting symmetry breaking 
inside the (a-GeTe(lll) triple-layers is responsible for the bulk Rashba-splitting. Of importance 
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for the surface-sensitive UV-SARPES on unpoled samples is that the spontaneous polarization 
of FE domains persist in the sub-surface region, as this ensures a net spin polarization also for 
the bulk states (SI IL3). 

Due to the small probing depth of ARPES, measured energy-, angle- and spin-resolved spec¬ 
tra can hardly be compared directly with a calculated band-structure because photoemission 
intensities do not relate directly to dispersions of bare bands. The task to describe the surface 
electronic structure differs from that in the bulk in several ways. Due to the broken crystal sym¬ 
metry the bulk spin-resolved band-structure calculations is different compared to a semi-infinite 
surface [Fig. l(f-g)], mainly due to the presence of surface states (SS). The star-like feature in a 
constant energy map (CEM) from bulk initial states orients along PK, whereas surface-like ini¬ 
tial states orient along PM. Furthermore, for certain binding energies and momenta, the canted 
spin of the bulk states is opposite to what is predicted for surface states with rather tangential 
helicity. Thus, mixing the spectral weights from surface and bulk states inevitably result in a 
competition of spin textures between bulk and surface states which in our case facilitate their 
discrimination in a spin-resolved experiment. 

Since ARPES data critically depend on photoemission final states, to build a bridge between 
SARPES and theoretical band-structure calculations we use so-called one-step photoemission 
model (ISM, see Methods). In this most rigorous interpretation besides bulk states also surface 
evanescent (SS) and surface-bulk-resonance (SBR) states are considered. To better understand 
their spectral weights, we employ the ISM determinant criterion nans] in Fig. l(i) in order 
to qualitatively discern their contributions in the measured valence states. Within the deter¬ 
minant logarithmic scale small values correspond to surface states with the maximum of their 
wavefunctions between the first atomic layer and the surface potential. States which embed 
in the bulk continuum have values close to one and are shown by black shades. On the other 
hand the SBR states, being partially reflected from bulk states, are given as red shades. The 
SS/SBR and SBR/bulk transitions appear gradually, but as a whole the SBR is the primary 
spectral weight dispersing across the gap, forming metallic states at Ep, and degenerating with 
bulk Rashba bands at lower binding energies (Ep). 

We first concentrate on the localization of the Z-point in k-space. The relevant region of 
fcz-momenta manifesting the Rashba-splitting is expected along the ZPZ direction, denoted in 
red in Fig. 1(a). The corresponding band maps shown in Fig. 2(a) have pronounced spectral 
intensity of the surface states at Ep and the Dirac point (DP), which we used to relate exper¬ 
imental data to band-structure calculations. The DP has minimum of -0.25 eV Ep around a 
hu of 22 eV, constituting the Z-point along the ZPZ direction. The VBM at the Z-point is ?^-70 
meV Eb in MPM and ?^-50 meV Ep in KPK mirror plane, as seen in Fig. 2(a,f). 

What is the origin of the metallic states at Fp? We note that SBR materialize as a super¬ 
position of bulk and surface spectral weight contributions, readily identified in corresponding 
theoretical OEMs in Fig. l(f-g). We note that by mixing SS with bulk states the SS typically 
lose their well-defined surface character at Fp by forming a SBR. Our results suggest that 
hybridization between the SS and bulk states is in-fact ’’mediated” by SBR. Recently such hy¬ 
bridization mechanism was also confirmed in Sb 2 Te 3 topological insulators m- An important 
observation is that the SBR fill the gap from the valence band maximum (VBM) all the way up 
to Fp, with progressively increased bulk-like spectral weight toward the VBM [Fig. 2(b)]. The 
in-gap states reaching Fp display rather 2D than 3D character in /cz, as seen in the band-map 
section I in Fig. 2(a). The metallic states filling the gap couple more to surface than bulk states 
for two reasons. First, consistently with the ISM predictions in Fig. l(i), the spectral weight 
at Fp for momenta fcp <0.3 A belongs to the SBR states and not to pure surface states. 
This corroborates their momentum distribution dispersing toward the pure surface states near 
fcp ?^0.4 A~ , as seen by the momentum distribution curves (MDC) in Fig. 2(e) and Fig. 3b. 
Second, the SBR states filling the gap clearly detach from the main 3D-valence band dispersive 
trend inside the band-map section II. 
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Figure 2: ARPES dispersion and comparison with semi-infinite surface calcula¬ 
tions. (a) UV-ARPES data taken along the ZFZ direction for selected photon energies point 
out the Z-point around 22 eV with the Dirac point (DP) located around -0.25 eV E^. The 
dashed lines highlight the dispersion of the DP along with the top of the valence band, indicat¬ 
ing a bulk gap of few tens meV. Yellow arrows indicate the 2D-like dispersion of the SBR inside 
the gap (section I), in contrast to the 3D-like dispersion of the bulk-like band in section 11. (b) 
OEMs for selected inside the gap measured at hz/ of 22 eV and 475 eV in (c). (d) ARPES 
ISM calculation calculated for hz/ of 22 eV along MPM showing the influence of matrix-element 
effects on the spectral weight, which expose the surface states in the band-map section III. (e) 
MDCs showing the 2D-like dispersion of SBR states at Ep. (f) Comparison of the semi-infinite 
and 10-layer slab calculations and corresponding second-derivative ARPES data in panel (g) 
highlights the SBR buildup in the valence band [yellow dashed lines in (f) and white dots in 
(g)]- 

To further understand our ARPES data it is useful to compare them with the semi-infinite 
and finite slab of 10-layers. As seen in Fig. 2(f), it is evident that the SBR-states pile up as 
a superposition of bulk-projected states giving rise to enhanced spectral intensity. Such SBR 
spectral intensity buildup is highlighted with dashed lines and is dominant in the valence band 
ARPES. It refers to SBR states with strong coupling to the bulk continuum. Consequently, in 
agreement with the ISM calculations, whole portions of pure bulk-VB states, denoted with red 
and blue lines, are overlaid in ARPES with the SBR states making their identification difficult 
by means of ground state bulk band-structure calculations mm- On the other hand the com¬ 
parison between the ISM ARPES calculation for hz/=22 eV in Fig. 2(d) with the corresponding 
experimental data in panel (a) is rather convincing. 

We discuss next our spin-resolved measurements near Ep. Spin-resolved populations of 
spin-up and spin-down electrons were measured along the r direction perpendicular to the scat¬ 
tering plane [Fig. 1(e)]. We use a well-established quantitative 3D-vectorial analysis [18] for 
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Figure 3: Energy dispersions measured by ARPES and SARPES near Ep. (a,b) (top) 
Valence band ARPES data with populations of spin-up and spin-down electrons near Ep, as 
a function of electron momentum along KFK and MrM(mid panels) [ with corresponding spin 
polarizations Px and Py at Ep (bottom panels), (c) Expected spin-resolved Eermi surface map 
(ESM) compared with experimental data in (d). Spin-resolved theoretical predictions from 
semi-infinite surface in Eig. 1(h) (gray arrows), are overlaid with smoothed x,y data from a,b 
(light-green arrows). Dark-green arrows mark the results from the 3D vectorial spin analysis 
based on simultaneously fitted MDC peaks [thin lines in (d), see SI V]. The difference in spin- 
texture between the simplified model sketched in panel (c) and the experimental result is due 
to spin hybridization of the SBR wavefunctions indicated by dashed arrows in panel (e) (see 
text). 


simultaneous peak fitting of spin-integrated MDCs and spin polarizations. The polarization 
curves are modeled until the best fit is reached by simultaneously fitting the MDC intensity and 
the polarizations Py and Pz (SI III). Of importance are the MDC-peaks 1-2-3 labeled in 
Eig. 3(a,b), the visibility of which is comprehensively described within the simultaneous fitting. 
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The resulting spin-texture at Ep is seen in Fig. 3(c). Thick arrows compare the 3D-fits with 
the in-plane polarization predictions x,y) from the semi-infinite surface (quiver plot in gray 
color). As a whole the SARPES data at Ey suggest that the SBR metallic states filling the 
gap have unconventional spin texture with parallel spin orientations different from the simple 
Rashba model. Our conjecture is that such spin topology in the momentum distribution can be 
reconciled via SOC-induced hybridization of the SBR-states. Highlighted in orange are regions 
where the bulk and surface contributions to the SBR wavefunctions hybridize. The resulting 
gap opening in momentum distribution, sketched in Fig. 3(e), is inducing a small spectral weight 
shift which realize the cochiral spin texture (SI III). Similar interband gap opening and spin 
reorientation has been observed in 2D Rashba systems [niEn]. As a matter of fact, the in-gap 
states inside the band-map section I in Fig. 2(d) display rather 2D than 3D character. 

To further examine our ARPES data we note that our experimental geometry predicts in 
ISM calculations a suppression of spectral weight in MFM below the DP for positive momenta. 
In fact, on the experimental side the bulk contributions to the SBR function wash out; un¬ 
covering the surface states inside the non-dispersive band-map section III [Fig.2 (a)]. This 
photoemission-transition matrix-element effect provides a direct evidence that the bulk and 
surface states indeed overlap in the initial states. As seen in Figure 2, a side effect of the 
vanished spectral weight intensity is that the OEMs appear with pronounced 3-fold rather than 
6-fold symmetry [2TJ [TT] . 
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Figure 4: Rashba splitting and out-of-plane spin polarization of bulk bands, (a) 

ARPES constant energy map (CEM) measured in the Z-point {hu =475 eV) at -0.4 eV Eb- 
Spin-resolved MDCs are measured along MFM (vertical) and KFK (horizontal) measured in 
equivalent Z-point at hv =22 eV. The spectral weight of bulk Rashba states for positive momenta 
vanishes as seen in Fig. 3(a), exposing the surface states (SS) denoted with violet arrows. The 
SS helicity (violet arrows) is opposite to pure Rashba bulk states (arrows with red and blue 
hues). Yellow frames highlight momenta where the ISM predicts pure bulk states (yellow 
arrows in Fig. l(f,i). The measured in-plane x,y polarization for both bulk bands aligns with 
canted spins consistently with the bulk Rashba initial states, showing a giant spin-splitting 
?^0.11 A . Along KFK the semi-infinite surface quiver is displayed in gray arrows. It 
is compared to measured (light-green quiver) and 3D-fitted (dark-green) spin textures, (b) 
Hexagonal warping in the out-of-plane spin polarization Pz measured along the paths (A) and 
(B). (inset) Corresponding ISM Pz calculation. 
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We now look at the influence of the opposite helicity of SS and SBR states on the mea¬ 
sured spin texture. It turns out that an ideal locus of k-space momenta for measuring pure 
bulk-Rashba properties is around -0.5 eV Eb in MFM . The black shades in the photoemission 
determinant along the green arrow in Fig. l(i) predict bulk states with distinct canted helicity 
for in-plane x,y polarization as seen in the inset in Fig. 1(f)]. The same canted 

helicity is perfectly reflected in the measured polarization in Fig. 4(a). Such experimental con¬ 
firmation traces back the magnitude and direction of the spin polarization vector to the states 
from which these electrons originate. This is non-trivial, since especially for systems with large 
SOC, the total angular momentum is relevant for photoemission process rather than only the 
spin quantum number [22]. Since the degree of spin polarization can depend strongly on the 
hz/ 122 ], the same states were probed with different photon energies (SI III). Because the same 
canted spin arrangement was confirmed in all cases, we conclude that the measured spin signals 
trace back to initial state bulk Rashba spinors, making this the first experimental confirmation 
of the spin texture of a FERS. 

Finally we comment on the out-of-plane polarization. The band anisotropy between FM and 
FK induce rather strong hexagonal warping effect responsible for the out-of-plane Pz spin po¬ 
larization component [TlESlEj. SARPES data in data seen in Fig. 4(b) confirm the hexagonal 
warping measured in two equivalent MFM directions denoted (A) and (B). As seen in the inset 
ISM calculations, the measured warping is in excellent agreement with measured data. 

In a broader perspective, the unveiled (a-GeTe(lll) bulk spin structure has far reaching 
consequences both for fundamental physics and applications. Based on our results (a-GeTe(lll) 
is the prime candidate in spin-orbit-driven Rashba physics in 3D k-space. Resolving the spin 
structure at Fp and inside the narrow gap is relevant in terms of potential device-oriented ap¬ 
plications. Namely, GeTe appears to have always p-type metallic transport properties due to 
Ge-vacancies which lead naturally to a large density of holes in the valence band states [25|. 
Nevertheless it manifests macroscopic, robust and reversible ferroelectricity promising for tech¬ 
nological applications. Our findings confirm that a bulk Rashba effect and ferroelectricity coexist 
and are coupled. We emphasize that the giant spin-splitting AkR ps 0.11 is formed by bands 
which interconnect the Rashba-split conduction- and valence band continuum. We anticipate 
that the (a-GeTe(lll) Rashba physics is likely to extend to Gei_^Mn^Te multiferroic Rashba 
semiconductors because the rhombohedral distortion persists up to x = 0.25 [Hj. Thus the 
quest for semiconductors with combined functional properties might extend to systems where 
magnetism, ferroelectricity and Rashba physics would offer unique property combinations for 
potential device applications. 

Methods 

Sample preparation method. 

Experiments were performed on 200 nm thick o-GeTe(lll) films grown by molecular beam epitaxy 
on BaF 2 (lll) substrates [121 US). A protective amorphous Te-cap used to avoid surface oxidation and 
degradation was removed in the ultrahigh vacuum chamber by annealing the thin film for 30 min at 
250 °C. The surface quality checked by FEED shows the 1x1 diffraction pattern of the o-GeTe(lll) 
surface seen in Fig. Ic. 

ARPES and spin-resolved ARPES 

The VUV-ARPES/SARPES experiments were performed at the GOPHEE end-station of the SIS beam¬ 
line [26], Swiss Light Source (SLS), using p-polarized photons in the energy range 20-25 eV, and an 
Omicron EA 125 hemispherical energy analyzer equipped with two orthogonally mounted classic Mott 
detectors. The whole set-up allows the simultaneous measurements of all three spatial components of the 
spin-polarization vector for each point of the band structure. The geometry of the experimental set-up is 
schematically shown in Fig. le. The spin-resolved MDGs were measured by rotating the sample azimuth 
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(j) such that MFM or KFK were aligned perpendicular to the scattering plane. Direct and immediate 
visualization of the spin splitting as a function of electron momentum was measured by exploiting the 
manipulator tilt angular degrees of freedom [see r in Fig. 1(e)]. The spin-integrated Fermi surfaces were 
measured with the channeltron detectors of the SARPES set-up by on-the-fly sweeping the manipulator 
tilt angle (r) for selected polar angles 0 as seen in the experimental setup in Fig. 1(e). The angular and 
resolution-broadened energy range in spin-resolved mode were 1.5°and 60 meV, respectively. In spin- 
integrated mode the resolutions were set to 0.5°and 20 meV. To enhance the bulk sensitivity compared 
to VUV-ARPES, experiments were also performed in the soft-X-ray energy range at the ADRESS beam- 
line m- The experimental results in the main text were reproduced on three samples with individual 
annealing preparation over different experimental runs. All data were taken at 20 K. 

First principle calculations 

The ab-initio calculations are based on density functional theory (DET) as implemented within the mul¬ 
tiple scattering theory (SPRKKR) [28l|29]. SOC has been naturally included by use of a fully relativistic 
four-component scheme. As s first step of our investigations we performed self-consistent calculations 
(SCE) for 3D bulk as well as 2D semi-infinite surface of (a-GeTe(lll) within the screened KKR for¬ 
malism [29] . The corresponding ground state band structures are presented in terms of Bloch spectral 
functions (BSE). The converged SCE potentials served as an input quantities for our spectroscopic in¬ 
vestigations. The ARPES calculations were performed in the framework of the fully relativistic one-step 
model of photoemission [ig ED] in its spin-density matrix formulation m, which accounts properly 
for the complete spin-polarization vector, in particular for Rashba systems like GeTe. Together with a 
realistic model for the surface barrier potential, one-step calculations based on a semi-infinite half-space 
configuration were decisive to substantiate the (a-GeTe(lll) surface-generated spectral features on both 
qualitative and quantitative levels. 
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